Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 19 January 2013 (MN WQi style file v2.2) 



Cosmography with strong lensing of LISA gravitational wave 
sources 



o 

Oh 
< 



o 

u 

43 
6 

o3 



> 

On 

o 



IS 



M. Sereno 1 ' 2 *, Ph. Jetzer 3 , A. Sesana 4 , M. Volonteri 5 

1 Dipartimento di Fisica, Politecnico di Torino, Corso Duca degli Abruzzi 24, 10129 Torino, Italia 
2 INFN, Sezione di Torino, Via Pietro Giuria 1, 10125, Torino, Italia 

3 'lnstitut fiir Theoretische Physik, Universitdt Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland 

4 Max Planck Institute for Gravitationalphysik (Albert Einstein Institute), Am Muhlenberg, 14476, Gohn, Germany 

5 Department of Astronomy, University of Michigan, Ann Arbor, MI 48109, USA 



19 January 2013 



ABSTRACT 

LISA might detect gravitational waves from m ergers of massive black hole binaries strongly 
lensed by intervening galaxies (ISereno et al.l 120 1 Qb . The detection of multiple gravitational 
lensing events would provide a new tool for cosmography. Constraints on cosmological pa- 
rameters could be placed by exploiting either lensing statistics of strongly lensed sources or 
time delay measurements of lensed gravitational wave signals. These lensing methods do not 
need the measurement of the redshifts of the sources and the identification of their electro- 
magnetic counterparts. They would extend cosmological probes to redshift z < 10 and are 
then complementary to other lower or higher redshift tests, such as type la supernovae or 
cosmic microwave background. The accuracy of lensing tests strongly depends on the for- 
mation history of the merging binaries, and the related number of total detectable multiple 
images. Lensing amplification might also help to find the host galaxies. Any measurement 
of the source redshifts would allow to exploit the distance-redshift test in combination with 
lensing methods. Time-delay analyses might measure the Hubble parameter Hq with accu- 
racy of > 10 km s _1 Mpc _1 . With prior knowledge of Hq, lensing statistics and time delays 
might constrain the dark matter density {5£Im ^ 0.08, due to parameter degeneracy). Inclu- 
sion of our methods with other available orthogonal techniques might significantly reduce the 
uncertainty contours for Om and the dark energy equation of state. 
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1 INTRODUCTION 

Observation of gravitational waves (GWs) by extragalactic sources 
is going to open a new window for as tronomy. The space-based 
Laser Interferometer Space Antenna (IDanzmann & et all Il996l 
LISA) is expected to observe up to s everal hundreds of events per 
year dSesana et al.ll2005ll200ll2010h . The loudest signals at LISA 
frequencies, / ~ mHz, should originate from coalescing massive 
black hole binaries (MBH Bs) with total masses in the rang e 10 3 - 
10 7 Mq out to z ~ 10-15 dHughesll2002l ; lKlein et all20Qgf) . 

Whenever a new experimental set-up to observe the uni- 
verse starts working, new possibilities open out. In a previous 
paper, we discussed a potential new channel for LISA science: 
multiple imagin g of GW s ources by intervening strong lensing 
galaxies dSereno et al.l 2010). Lensing of distant so urces has been 
long considered as a test fo r cosmological t heories dRefsdall 
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erences therein). GW sources might allow for a variation of these 
classical tests. 

The main novelty of making cosmography w ith LISA rel ies on 
the property of MBHBs of being standard sirens (ISchutz 1986). The 
luminosity distance to the inspiral GWs can be determined with 
good accuracy and sev eral methods have alr e ady been proposed to 
explo i t this property ijHolz & Hughes 2005 : Jyan den Broeck et al] 
l20ld : IShapiro et alj|201Ct iHilbert :et aljkoich . The main idea on 
the table is to build-up the Hubble diagram. The relation between 
distance and redshift changes for different cosmological theories or 
different cosmological parameters. However, the redshift cannot be 
measured from the analysis of gravitational waves alone. The use 
of MBHBs as cosmological probes should rely on the identification 
of the electromagnetic counterpart in order to measure the redshift 
of the source. 

Lensing offers an alternative tool. In the classical Hubble dia- 
gram, we compare observed distances with theoretical expectations 
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depending on cosmology and redshift, which has to be measured 
in an independent way. When performing cosmography with lens- 
ing, we compare some quantity inferred from the lensing analysis, 
which is in general a combination of cosmological distances, with 
the theoretical expectation, which depends on the cosmological 
model. To perform these tests we need either the measured distance 
to the source or the redshift. In usual lensing studies of quasars 
or radio-sources, we exploit the measured redshifts of the sources 
l lKochanelj[T993l ; IChaell26o3) . In lensing studies of LISA sources, 
we should determine the distances to the sources. The identifica- 
tion of the electromagnetic counter-part of the signal is then not 
necessary to perform the cosmological tests and we can bypass the 
weakest link in the already proposed cosmographic methods with 
LISA. 

The chances of observing multiple images of the same GW 
source with LISA are sizeable for a broad variety of formation his- 
tories, ranging from < 40% to < 100% for a 5 years mission ac- 
cording to the redshift distribution of th e sources and their intrinsic 
signal to noise ratio l lSereno et al.luOlOl) . It is then worthy to inves- 
tigate what kind of cosmographic tests might be possible with LISA 
lensing. Strong lensing by ground-based GW detectors was dis- 
cussed in lWang et alj d 1 996 >. that considered the constraints on the 
amount of matter density in compact lenses as derived from lensing 
statistics in the context of advanced LIGO type detectors. Gravi- 
tational lensing of GWs was also considered to measur e the rela- 
tive tr ansverse velocity of a source-lens-observer system dltoh et all 
l2009h or to obtain inf ormation about the typical mass of lens ob- 
jects jYoo et al.ll2007» or as a systematic effe ct hamperi ng the de- 
tection of a very weak primordial GW signal (Seto 2009). 

Here, we focus on the LISA mission and its potential to probe 
cosmological parameters by strong lensing methods. In Sec. [2] we 
review some basics of lensing of GWs by galaxies. In Sec. [3] we 
develop a formalism for the computation of probabilities to detect 
transient lensing events. In Sec. [4] we discuss how to infer the lu- 
minosity distance to a multiply imaged GW source. Cosmological 
tests based on either lensing statistics or time delay measurements 
of GW sources are introduced in Sec. [5] Section|6]lists the assump- 
tions made for a plausible lensing scenario whereas Secs.|7]and [8] 
contain our forecasting of the lensing test accuracy. Section [9] is 
devoted to some final considerations. 

As reference model we consider a flat ACDM model with 
f^M = 0.3, Q.a = 0.7, and h = 0.7 where h is the Hubble constant 
H in units of 100 km s _1 Mpc -1 . 



2 BASICS OF GRAVITATIONAL LENSING 

The statistics of gra vitational lenses have become a standard 
tool for cosmology dKochaneld [19961 ; IChad [20031 ; ISerenol 120051 ; 
IZhu & Serenol2008[ and references therein). The differential prob- 
ability of a background source to be lensed by a foreground galaxy 
with velocity dispersion between a and a + da in the redshift in- 
terval from Zd to Zd + dzd is 



d 2 t dn . . , . cdt . . 
-i — -7- = -r(Zd,T)s CI (Zd,a) — (Zd), 
dzdda da dZd 



(1) 



where s cr is the cross section of the deflector and dn/da is the 
differential number density of the lens population. 



2.1 Lens mass density 

Departures from spherical symmetry and details of the radial mass 
distribution of the lens galaxy induce a relatively small effect on 
len s statistics and are unimportant in altering the cosmologica l lim- 
its dMaoz & Rixll 1991 iKochaneiill 19961 ; [Mitchell et alj|20bl) . We 
can therefore approximate early-type galaxies as singular isother- 
mal spheres (SISs). Two images of a compact source aligned with 
the lens form at x± — y ± 1 if y < 1, with flux magnifica- 
tion p± = (1/y) ± 1. Here, x and y are the image and the 
source angular position normalised to the angular Einstein radius, 
8 E = 4ir(a/c) 2 D dB /D E . D d , D dB and D E are the angular diameter 
distances between the observer and the deflector, the deflector and 
the source and the observer and the source, respectively. 

In GW lensing, at variance with the usual lensing, we directly 
observe the wave-form, which is amplified by A± = y/JI±. The 
delay between the arrival time of the images, At = i_ — t+, is 



At = At z y,At z 



15 days 



32tt 2 fa\ i D d D. 



r \c; D B 
D d D ds 



200km/s J D s x 0.2Gpc 



(1 + zd) (2) 
1 + z A )y. 



2.2 Lens population 

The lens distribution can be modelle d by a modified Schechter 
function of the form jSheth et al.l2 003) 



dn 
da 



(3) 



where a is the faint-end slope, f3 the high-velocity cut-off and n, 
and <7» are the characteristic number density and velocity disper- 
sion, respectively. An evolving galaxy density can be paramet erised 
with ruQz) = n», (l + z) 3- ""* and a, = cr* j0 (l + z)""* dChad 
120071) . For a constant comoving number density, v n * = v a * = 0. 



2.3 Detection thresholds and bias 

Lens dis covery rates are affected by the abil ity to observe multiple 
images dKochanelal 19931 ; ISereno et alj|20101) . For optically lumi- 
nous quasars or radio sources, lensed systems are detected by se- 
lecting resolved multiply images; intrinsic variations show up with 
a time delay in each image. On the other hand, due to the low angu- 
lar resolution and the transient nature of the sources, lensed GWs 
in the LISA context are detected as repeated events in nearly the 
same sky position. 

The source position is limited to an allowed range, y m i n ^ 
y ^ J/max, for which multiple images are detectable. The upper 
limit j/ max depends on the lens mass, the arrival time, the threshold 
signal to noise r atio (SNR t h) and t he unlensed amplitude of the 
source (SNR int ) dSereno et alfcoioh . 

We require that lensing amplification pushes the signal of the 
second image above threshold, A- > SNR^/SNRmt, which lim- 
its the source position to 



y ^ ?/max 



'SNR-/SNR in 



+ 1 



(4) 



The minimum y m i n excises the region near the central caustic 
where wave optics is effective and the interference pattern covers 
the multiple images. For the LISA wave-band, geomet ric optics is 
adequate and we can put y m in = dSerenoetal.l2010l) . 

For our forecasting we deal with known properties of the 
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source population so that we do not need to correct for the 
magnification bias {Fukugita & Turner 1991 ; Fukugit a et al .11 199*3 ; 
ISereno et al.l2010h . 



3 STATISTICS FOR TRANSIENT PHENOMENA 

Due to the finite duration of the survey, lensing statistics of transient 
pheno mena involve acc ounting for some missing events due to time 
delay (Ogu ri et al.ll2003h . The cross section of a SIS for a time- 
limited survey is 



4 DISTANCE DETERMINATION 



(2d, c) = t-Re 



f(At(y;a, z A ))ydy, 



(5) 



where _Re = D^Oe is the Einstein radius and /(At) is the frac- 
tion of lenses with time delay At that can be observed. For lasting 
images, as for quasars or radio-sources, /(At) = 1, 



2/min) 



(6) 



( 2 

I 2/ max 



J/min) 



I 3 
I l/max 



3 > 
l/min) 



(8) 



We consider a uniform distribution of arrival times during the 
time survey, T SUI . If the monitoring is continuous 

/(At) =l-p- (7) 

J sur 

for At < T sm . and otherwise (Oguri et al. 2003). The resulting 
cross section weighted for the arrival time distribution is then 

2 At z 

3 T sur 

In the following, we specialise to the case of geometric optics 

(j/min = 0). 

The differential probability dr / dz A for a source to be lensed 
by a deflector in the interval between Zd and Zd + dzd can be ob- 
tained integrating the differential probability in Eq. Q]. For each 
lens redshift Zd there is a maximum velocity dispersion o- max 
such that Scr(c) < for a > cr max . In practice, <r max is quite 
large (3> and the corresponding galaxy density is almost null, 
dn/ da {a max ) ~ 0. When integrating the differential optical depth, 
we can then take <r max — > oo. We get 

r \ 4 + a ] 

— ~ n < ., s CI .(a-,,o;zd) —. r l 1 - 1 -^'^"^ — 

dz d ii(zd) 



r [f] 



_ 2 At(o-« |0 ,y n 

3 isur 



'(l + 2d) 



(9) 



where Y is the Euler gamma function. 

The total optical depth for multiple imaging of a compact 
source, r, gets a quite compact form for a constant comoving den- 
sity, v n » = v a * = 0. After integration, 



T=^[D B (l + Zs)T „ 



2 

max 



where 

F* = 



167T n, o 



At. = 32 



c 

4 D, 



i r[(8 + a)//3] At* 
7r[(4 + a)//3] T sur 

4 T [(4 + 



r[a//3] 

vr I— ^— I — (1 + z s )y max . 



(10) 



(11) 



(12) 



The optical depth r is a function of the source redshift z s and in- 
trinsic SNR (through j/ max ); the cosmological parameters enter in 
the angular diameter distances. 



In a seminal paper, ISchutzl ( ll986l) showed that measurements of the 
amplitude, frequency and frequency derivative of the inspiralling 
massive binary black holes could yield a precise estimate of their 
luminosity distance. This would allow us to measure the distance 
in a novel way, making GW sources potentially powerful standard 
sirens. 

Each harmonic of the inspiral polarisations is inversely pro- 
portional to the distance, p oc 1/Dl, where Dl is the lumi- 
nosity distance. For well modelled systems, LISA will be able 
to measure the luminosity distance (but not the redshift) to mas- 
sive BH binaries with 1-10% accuracy. The measurement preci- 
sion is largely limited by pointing error and weak lensing distor- 
tion dLang & Hughesll2006l : iKlein et alJl200Sl) . Gravitational lens- 
ing will randomly magnify or demagnify MBHB signal, and thus 
system atically modify any distance measurement dHolz & H ughes 
l2005t) , 

LISA should detect independently two above threshold signals 
in the same sky position in order to claim lensing. The two lensed 
wave-forms are identical, apart from an overall factor connected 
to lensing amplification. For each image p± oc A±/Di,, where 
the suffix refers to the lensing parity of the image and not to the 
polarisation. 

The distance to the lensed source can be inferred in the follow- 
ing way. The ratio of the amplifications A- /A+ can be measured 
from p-/p+, and one can directly infer the source position, y, 



l + \A^/A, 



(13) 



Once the source position is known, one can quantify the amplifica- 
tion of each image, A± and estimate the luminous distance. As a 
consistency check of the lensing hypothesis one should find that the 
value of y does not depend on the frequency as is the case instead 
for the amplitudes. Since lensing amplification can be expressed in 
terms of the scaled lens position y, knowledge of the lens velocity 
dispersion is not needed. Such methodology can be easily gener- 
alised to lenses more complex than the simple SIS. 

The method above is effective only with regard to strong lens- 
ing by a single deflector plane. Weak lensing is well recognised 
as a potential noise in the determination of the distance to GW 
sources. Amplification due to large scale structure cannot be fil- 
tered out and will still contribute the main uncertainty in the de- 
termination of Dl- There are proposals on how this effect could 
be, at least partially, corrected. Convergence maps reconstructed 
using galaxy flexion in addition to shear mig ht help to reduce the 
lensing-induced di stance errors by up to 50% dShapiro et al.ll20ic| ; 
Hilbe rt et alJkoiCh . A Gaussian distance error with a standard de- 
viation of 10% was then added to each luminosity distance to simu- 
late the effect of weak lensing errors. A more realistic error should 
grow with redshift, but we also tested that results are nearly unaf- 
fected by assuming a 20% error. 

Magnification o r demagnific ation of signals affect strong lens- 
ing observables too dAsadal 1998). However, the universe is quite 
homoge neous at the high redshifts (z > 10) probed by LISA 
sources dWandll999h . Magnification effects are also somewhat 
washed out after av eraging over many independent lines of sights 
dSereno et al.l 120021) . Lensing statistics are then not heavily af- 
fected by the lensing dispersion in the distance-redshift relation 
dCovoneetai]|2005l) . 
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5 COSMOLOGICAL TESTS 

The lensing optical depth depends on both source redshift and dis- 
tance. In standard lensing statistics of radio sources or quasars, the 
source redshift is known and we can estimate the distance assum- 
ing a given cosmological model, i.e., a given set of cosmological 
parameters. Dealing with lensing statistics of observed GWs, we 
know the luminosity distance to the source and we have to work 
the other way: given a set of cosmological parameters, we estimate 
the source redshift and, on turn, the angular diameter distance. 

The analysis of the GW-form allow us to infer the luminos- 
ity distance to the source, whereas lensing quantities are written in 
terms of angular diameter distances. In the general theory of relativ- 
ity, luminosity and angular diameter distances are related through 
the Etherington principle, Dl = (1 + z) 2 D. Throughout, lumi- 
nosity distances are denoted by the suffix "L". It is also useful to 
remember that in a flat universe D& s = D s — Dd (1 + / '( 1 + z B ). 

5.1 Lensing statistics 

The likelihood function for lensing statistics ca n be written as 
( lKochanelj[T9^ : IChaell2003l : lMitchelle7aIll2005l) 

N V N h 

£ s tat(fii) = JJ(1 - n) Y[pi,h (14) 

where {f2, } is the set of cosmological parameters that characterises 
the model under i nvestigation. Len s statistics do not depend on the 
Hubble constant iKochanekl 1 19931 : IChad 120031) . ]Vl is the num- 
ber of multiple-imaged sources and iVu is the number of unlensed 
sources, pij is either dr/dz^ for known lens redshift or r for un- 
known Zd- In principle, a test based on lens statistics can then be 
performed even without lens identification. 

For given values of the cosmological parameters, the relation 
between z and Dl is unique. Given the measured Dl, we can then 
compute the angular diameter distances and £ s tat for each set of 
cosmological parameters. 

The main source of statistical uncertainty is by far the Pois so- 
man noise due to the small number of events. For our forecasting 
about the accuracy of lensing statistics, the error on Dl, which is 
of order of < 10%, is then negligible. 
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Figure 1. Probability density functions for the redshift of the 3 lensed 
sources. PDFs have been obtained by repeatedly extracting a triple from the 
intrinsic redshift distributions of the simulated LISA sources. Each source 
was weighted by its optical depth to lensing. 




Figure 2. Differential lensing probability for different source redshifts. 
Probabilities were obtained through Eq. l|2}. Intrinsic SNR are 400, 12 or 
36 for z s = 4.5, 6.2 or 9.3, respectively. 



5.2 Time delay 

It has been long known that time delay measurements c an be used 
to co nstrain cosmological parameters (Refsdal 1966; Sa ha et all 
2006). If we measure At and constrain the lens properties with 
either lensing or other follow-up imaging or spectroscopic obser- 
vations, we can get an estimate of Da* = -Dd-Dds/-Ds and in turn 
constrain the cosmological parameters. Differently from lens statis- 
tics, we need the lens redshift to carry out cosmography with time- 
delays. For each measured At, we can write a \ 2 contribution as 



D 



Obs 



Da* 



+ 



D 



SD Ob s 



(15) 



In the above equation, the apex "Obs" denotes the measured value 
of the corresponding quantity. D^ s is measured with the lensing 
analysis; the observed value of the luminosity distance D® hs is 
obtained from the wave-form investigation. D^t and Z?l are func- 
tions of the cosmological parameters and of the source redshift. 
Due to the uncertainty 5D® bs on Dp bs , the source redshift is not 
perfectly known and has to be considered as a model parameter. 



The second term in the right-hand side of Eq. ( 115) shapes the in- 
formation on the luminosity distance obtained from the analysis of 
the signal as a Gaussian prior. Once data analysis of LISA data is 
really performed, best estimates of luminosity distances to mergers 
are going to be determined together with their probability distribu- 
tions. More realistic priors for the time delay likelihood should be 
shaped after these computed distributions. 
The likelihood built on time delays is 



CAt(£li-,h;z s ,i) oc exp{-^ y^xlt,}, 



(16) 



where the sum extends on the lensed systems. The unknown source 
redshifts {a s ,i} are JVl additional model parameters. The cosmo- 
graphic approach with time-delays is not affected by Poissonian 
noise, whereas the error on Dl is an important source of statistical 
uncertainty. The main cosmological dependence is on the Hubble 
constant, since distances oc 1/Hq. 
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6 A SCENARIO FOR LENSING 

We tested the predictive power of the lensing methods on a class 
of flat cosmological models with dark matter, Qm, and a dark en- 
ergy fluid whose equation of state is parameterised by a constant 
equation of state ui. Within such context, our "true" ACDM model 
is described by Qm ~ 0.3, and oj — — 1. State-of-the art analy- 
sis of the cosmic microwave background radiation with informed 
priors can constrain t he Hubble constant with an accuracy < 2% 
( Kom atsu et ai] |201(t) and estimates should be even more tight by 
the time LISA is flying. Such error is much smaller than the uncer- 
tainties on Dg bs and _Dg^ s . 

We then consider two frameworks for cosmological parameter 
forecasting. Either we focus on the determination of the Hubble 
constant or we take Ho as known and concentrate on what lensing 
with LISA can do for the dark sector. 



6.1 Deflectors 

A proper modelling of the distribution of the lensing galaxies is 
central in lensing statistics. Early-type or late-type populations con- 
tribute to the lensing statistics in different ways and type-specific 
galaxy distributions are required. As a conservative approach, we 
did not consider lensing by spiral galaxies. Late-type galaxies con- 
tribute no more than 20-30% of the total lensing optical depth and 
the kno wledge of their number density i s plagued by large uncer- 
tainties dChaell2003MMitchell et all2005l) . 

In our analysis we used the results of Choi et al. (2007) who 
analysed data from the SDSS Data Release 5 to derive the ve- 
locity dispersion distribution function of early-type galaxies. They 
modelled the galaxy population as a modified Schechter func- 
tion with n»,o = 8.0xl0~ 3 /i 3 Mpc~ 3 , cr*, = 144±5 km s~\ 
a = 2.49 ± 0.10, and /? = 2.29 ± 0.07. We keep constant the 
comoving number density of galaxies. 

Even if the LISA angular resolution is quite poor, ~ 30'- 
1 deg, lenses are expected to be very massive and luminous. In 
case the deflector can be identified and its redshift measured, the 
corresponding Einstein radius can be estimated and one can accu- 
rately determine the angular position of the source. This way one 
could identify the source (more precisely its images) as well with 
follow-up observations and thus get also its redshift. 



6.2 Sources and lenses 

The number of detec table multiple ima ges depends on the build- 
up formation history dSereno et al.ll2010l) . Massive mergers at high 
redshifts can produce very loud GW emission with a noticeable 
optical depth to lensing but are expected to be quite rare. On the 
other hand, minor mergers might be more frequent but with a lesser 
lensing probability due to intermediate redshift and lower intrin- 
sic signal-to-noise. The total lensing probability balances the to- 
tal number of events, from f ew dozens to seve ral hundreds de- 
tectable coalescences per year dSesana et al J2010h . and their intrin- 
sic loudness. Throughout, the monitoring time for LISA is fixed to 
Tsur = 5 ye ars. The expected number of lensing events goes from 
< 1 to < 4 dSereno et aljioloh . 

We considered two scenarios for the abundance of LISA 
sources. In a first pessimistic case, we derived cosmological con- 
straints from a single lensing detection. In a more optimistic sce- 
nario, three multiple image events might be seen by LISA. This is 
the case for example for some hybrid formation histories, such as 



the "Hybridll" model of ISereno et al.l d2010h . Mergers were gen- 
erated in the reference ACDM model. For each source we com- 
puted the lensing probability. As a threshold for detection we took 
SNRth = 8. A mean of ~ 2.8 detectable lensing events is ex- 
pected. 

We extracted the properties of the three lensed sources from 
their parent distribution, weighting each source by its optical depth 
to lensing. The probability density function for the source redshifts 
(sorted according to their z) is plotted in Fig. Q] Accordingly to 
the mean properties of the distributions, we fixed the three source 
redshifts for our forecasting at z s ~ 4.5, 6.2 and 9.3. 

The corresponding lens redshifts were extracted as the Zd 
which maximises dr/dzd for a given z s , see Fig. [2] We got Zd — 
0.9, 1.0 and 1.0. Since the shrinking of the cosmological volumes 
with redshifts, the peak of the probability is always at Zd ~ 1, with 
no regard to the source redshift. Together with the expectation of 
being very massive, their moderate redshift might make the lenses 
identifiable. 

For the pessimistic scenario, we computed the most likely 
lensing configuration in the case of only one lensing event detected 
in the "Hybridll" formation history. The source is at z s ~ 7.4 and 
the corresponding deflector is at Zd — 1.0. 



7 FORECASTING FOR UNKNOWN SOURCE 
REDSHIFTS 

For our forecasting, we exploit different priors according to 
whether we focus on either Ho or Qm and u). When we deal 
with the Hubble constant, the accelerated expansion of the Uni- 
verse is assumed to be propelled by a cosmological constant, so 
that uj — —1 is our delta prior for the dark energy. As far as 
the dark matter is concerned, as a first case we assume no previ- 
ous information, i.e., a prior in the form of a uniform distribution, 
P(Qm) = 1 for ^ Qm ^ 1. Alternatively, we shape results from 
other methods as a very mild Gaussian prior centred on the "true" 
value Qm = 0.3 and with a quite large dispersion o~n M = 0.1. The 
prior on h is a flat distribution non null between and 2. 

When we focus on the dark sector, Ho is kept to its reference 
value. As prior for the cosmological parameters we consider a uni- 
form distribution in the square Qm ^ 1 and — 2 ^ uj ^ 0. 

7.1 Lensing statistics 

Lensing statistics of GW sources might constrain the dark matter 
density and the dark energy equation of state. It does not depend on 
Ho- We simulated the measurements of the luminosity distances by 
adding a Gaussian noise of 10% to the true value. However, since 
the Poissonian error is the main one, even by doubling the error on 
the distances the inferred posterior probability of Qm and uj from 
lensing statistics is unaffected in any sensible way. 

The likelihood was computed in a grid. The posterior proba- 
bility function for the cosmological parameters is plotted in Fig. [3] 
Lensing statistics can be performed considering either known or 
unknown lens redshifts. Precise information on Zd makes the upper 
limit on the estimated Qm much stronger. The accuracy on Om is 
~ 0.24 for 3 lenses with measured redshift. 

7.2 Time-delays 

The distance combination Da* is the basic source of information 
on cosmology from time-delay measurements. The uncertainty in 
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Figure 3. Posterior probability for Qm and u as derived from lensing 
statistics in the case of 3 lensed sources. Grey-shadowed regions or thick 
contours (full, dashed or long-dashed lines) refer to lensing statistics ei- 
ther exploiting information on lens redshift or assuming to be unknown. 
Contours are plotted at fraction values exp(— 2.3/2), exp(— 6.17/2), and 
exp(— 11.8/2) of the maximum, which would denote confidence limit re- 
gion of 1 -, 2- and 3-<r in a maximum likelihood investigation, respectively. 
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Table 1. Predicted uncertainty on the Hubble constant (column 5) under dif- 
ferent hypotheses for number of observed lenses (first column), uncertainty 
on D At (second column), previous knowledge of the lens redshifts (col. 3) 
and prior on Sljvi (col. 4; a Gaussian prior centered on 0.3 with dispersion 
of 0.1 or a uniform distribution between and 1). 



the measurement of D A t depends on the accuracy with which we 
can measure the time-delay and the lens mass distribution. For 
transient events, t ime delay can be measured with great accuracy 
(Oguri et al. 2003) so that the main source of error is from the un- 
certainties in the mass distribution of the deflector. Lens are ex- 
pected to be at redshifts Zd ~ 1 and to be very luminous so that 
follow-up spectroscopy and photometry should be performable. 
Further constraints come from the lensing modelling exploiting im- 
age positions and time-delay. An uncertainty on D&t of order of 
20% is a conservative assumption; an uncertainty of 10% is more 
optimistic and more adequate to the case of a deep follow-up. 

When considering time-delay, together with the cosmologi- 
cal parameters we have to consider as unknown model parameters 
the redshifts of the lensed sources. For our optimistic (pessimistic) 
lensing scenario, we have then to add 3(1) unknown source red- 
shifts to the cosmological parameters to characterise the model. 

7.2.1 H 

Forecasting for the Hubble constant is summarised in Table [T] 
The quoted uncertainty on Ho is the standard deviation of the 
marginalised posterior probability density. For each case, the pa- 
rameter space (Ho and Qm plus the source redshifts) was ex- 
plored running 4 Markov chains of 2.5 x 10 4 samples each. Chain 
convergence was checked through the Gel man and Rubin ratio 
( Gelm an & Rubirtl992tlLewis & Bridld2002h . which was well un- 
der 1.05 for each parameter. 

Even a mild prior on Q,m is very helpful in breaking the degen- 
eracy between Ho and the density parameters. With 3 lenses mod- 
elled with deep follow-up observations, an accuracy of Sh ~ 0.12 
can be achieved. This figure is larger than present uncertainty but 
it would provide a direct test on Ho without need for calibration or 




Figure 4. D& t (zd = 1,2s = 7) (thick lines) and D a (z s = 1) (thin 
lines) as a function of the cosmological parameters. Values of -Da* (D b ) 
are normalised to their value at = 0.3 and U3 = —1 and run from 0.7 
(0.8) to 1.6 (1.4) in steps of 0.1. The Hubble constant is fixed to h = 0.7. 



distance scale ladder in an unexplored redshift range, z < 10. Even 
a single lens detection could provide interesting limits, 8h ~ 0.20 
or 0.26 according to the a-priori degree of knowledge on the dark 
matter density. 
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Figure 5. Marginalised posterior probability for SIm and u as derived from 
3 measured time-delays. Grey-shadowed regions are for 5D9]? s /d9]? s = 
10%; thick (full, dashed or long-dashed) lines are for SD^/D^ 3 = 
20%. Contour values are as in Fig. [3] 



7.2.2 Dark energy 

Time-delay measurements might also constrain dark matter and 
dark energy. In Fig. H] Da* is plotted as a function of the cosmo- 
logical density parameters. Lensing by LISA exploits very distant 
sources. The dependence on cosmological parameters is then some- 
what orthogonal to tests exploring a lower redshift range, such as 
observations of type la supernovae. This can be seen by compar- 
ing constant contours of either Da* or luminosity distance in the 
region of interest of the parameter space. 

Since the total number of lenses is small, lensing methods by 
their own sample only a restricted redshift interval and are not able 
to break the degeneracy in the Qm-u plane. As can be seen in 
Fig. [4] even with an extremely accurate but single measurement 
of Da«, the degeneracy on the dark energy equation of state would 
be severe. A better observational accuracy could only reduce the 
confidence regions which would still be very elongated along the 
degeneracy direction. Even with a very small error on Da*, the ac- 
curacy on Qm is still affected by the degeneracy. The degeneracy 
can be broken only by combining with other tests exploring differ- 
ent redshifts. 

Together with Qm and uj we have to consider as unknown 
model parameters the redshifts of the lensed sources. For our opti- 
mistic lensing scenario, we have then to add 3 unknown lens red- 
shifts to the 2 cosmological parameters. We explored the parameter 
space by running four Markov chains of 5 x 10 4 samples each. 
This was more t han enough to reach convergence, with the Gel man 
and Rubin ratio dGelman & Rubinlll992l ; ILewis & Bridlel2002h be- 
ing well below 1.01 for each parameter. The prior for each source 
redshift was uniform and non zero for z s > Zd ■ 

The marginalised posterior probability function for Om and u 
is plotted in Fig. [5] After marginalization, the uncertainty on Qm 
is ~ 0.11 or ~ 0.21 for SD^/D^ = 10% or 20%, respec- 
tively. Since the redshift range is the same and the time-delay test 
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Figure 6. Marginalised posterior probability for S!m and ui as derived from 
both lensing statistics and time-delays in the case of 3 lensing events. Re- 
gions and contours are as in Fig. [5] 




Figure 7. Marginalised posterior probability for Um as derived from both 
lensing statistics and time-delays in the case of three lensing events. The 
thick and the thin lines are for <5r>°ks/£)Ob s = io% or 20%, respectively. 

is still based on angular diameter distances, the degeneracy in the 
parameter space is similar to lensing statistics. 

A less optimistic expectation for the accuracy in the determi- 
nation of the luminosity distances translate into worse cosmolog- 
ical constraints. However, the effect is not so dramatic. Assuming 
a 20% relative error on D^, the uncertainty on Qm after marginal- 
ization would be ~ 0.12 or ~ 0.22 for SD^/D'ffi* = 10% or 
20%. 

The time-delay test might be effective even for a sin- 
gle detected lensing event, when we get SQm ~ 0.18 for 

SD<£ S /D<£ S = 10%. 

7.3 Lensing statistics and time-delays 

To consider lensing statistics and time-delays at the same time, we 
used a combined likelihood C oc £ a tat x £a*- The parameter space 
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Figure 8. Posterior probability for Qm and u as derived from 3 measured 
time-delays of sources with known redshifts. Grey-shadowed regions are 
for 5DgJ?s/£)Ob s _ 1Q% . mick lines ^ for 5D Obs/£,Obs = 2 o%. 

Contour values are as in Fig. [5] 

was explored as in the time-delay case. The marginalised posterior 
probability function for Om and u is plotted in Fig.[6] We assumed 
that the Hubble parameter is known. The confidence regions are 
more confined than for lensing statistics or time delay alone, but the 
equation of state is still undetermined. The marginalised probability 
for JIM is plotted in Fig. [7] 

The dark matter density parameter might be determined with 
an accuracy of ~ 0.10 (0.18) for 3 lenses and SD^/D^ = 
10% (20%). Due to the small number of lenses, the final accuracy 
is mainly determined by the time-delay test. 



8 FORECASTING FOR KNOWN SOURCE REDSHIFTS 

In a best-case scenario source redshifts might be determined. In 
fact, lensing might help in finding the electromagnetic counterparts 
of the GW signals. Together with the GWs from the coalescence, 
galaxies harbouring the merging black holes should be lensed as 
well. Multiple optical images of the host galaxy might be then de- 
tected in the nearby of the deflector. The fluxes of the images should 
be magnified in agreement with the relative amplification of the 
GW signals, and their position relative to the lens should match the 
prediction from the analysis of the lensed wave-forms, see Sec. [4] 
Follow-up observations should make the photometric or spectro- 
scopic determination of the source redshift possible. In principle 
and very speculatively, the source redshift might be estimated even 
if only the most magnified optical or radio image of the host galaxy 
is brought over-threshold. In such a case, the lensed galaxies could 
not be found as usual by comparing the spectra or the multi-band 
fluxes of near galaxies, but we should look for a galaxy whose dis- 
tance from the lens is in agreement with the prediction from the 
lensing analysis of the lensed GWs. 

However, the inclusion of magnificati on ratios in lensing anal- 
yses is problematic for two main reasons dSaha et al1l2006l) : first, 
optical flux ratios may be contaminated by microlensing and dif- 
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Figure 9. Posterior probability for Qm and ui as derived from 1 measured 
time-delay of a source with known redshift. Regions and contours are as in 
Fig.S 
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Figure 10. Posterior probability for Cm and u) as derived from both lensing 
statistics and time-delays in the case of three lensing events with known 
source redshifts. Regions and contours are as in Fig. [8] 

ferential extinction; second, relative magnifications along different 
directions are weakly coupled with time delays, because magnifica- 
tion measures the local second derivative of the arrival time. Future 
searches for host galaxies should account for this. 

Once the redshift of a lensed source is known, one can 
strengthen the lensing tests by exploiting the distance-redshift 
comparison for the source, as usually done when building the 
Hubble diagram. The additional information on the source red- 
shift makes the lensing constraints on the cosmological parameters 
much tighter. To estimate the impact of known source redshifts on 
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our forecasting, we re-made the analysis of Sec.[7]by fixing the z a 's 
to their actual values. Now, the only free parameters to be inferred 
are the cosmological parameters. In our studio-case, we are then 
left with just two parameters (either Ho and JIm or Qm and oj). 

Constraints from the analysis of time delays get much 
stronger. The second term in the right-hand side of Eq. ( 1151 , which 
before stood to account for an uncertainty and could only help to 
estimate the source redshift, provides now an additional strong con- 
straint on the cosmological parameters. Posterior probability den- 
sities were computed on a two-dimensional grid. 

Uncertainties on h range from 0.1 to 0.3, see Table [T] Once 
the source redshift is fixed, even a single lensing event can provide 
very interesting bounds. The expected accuracy on Ho from lensed 
GWs compare s very well with more standard results from time- 
delay analyses. Saha et al. ( 2006) obtained 5/i ~ 0.1 for 10 lensing 
galaxies. 

Results for time-delay constraints on dark matter and dark 
energy are represented in Fig. [8] for the case of 3 lenses and in 
Fig. [9] for 1 lens. Thanks to the additional constraint from the 
distance-redshift relation, confidence regions in the Qm-oj plane 
are much more restricted than for unknown source redshifts. Due 
to the degeneracy, the accuracy improvement for the determination 
of Q,m is not so significant. With three lenses, the uncertainty on 
n M is ~ 0.08 or ~ 0.09 for SD^/D^ = 10% or 20%, re- 
spectively. Improvement is more significant for the case of only 
one lens, when the uncertainty on Qm is ~ 0.11 or ~ 0.12 for 
6D^ a /D^ s = 10% or 20%, respectively. 

The lensing statistics test is not very helpful in such an opti- 
mistic scenario. Combined posterior probability densities are plot- 
ted in Fig.[l0] Confidence regions are very similar to those obtained 
exploiting only the time delays. 



9 CONCLUSIONS 

The serendipitous discovery of multiple images of gravitational 
wave sources by LISA might offer new possibilities for astronom- 
ical investigations. Gravitational lensing has been long considered 
as a tool for cosmography. Here, we considered the peculiarities of 
cosmological tests with lensed GW sources. Both time-delay mea- 
surements and lensing statistics were investigated, which required 
to develop a treatment of lensing probabilities for transient events. 

A couple of features of the discussed methods deserves partic- 
ular attention. The main appeal of lensing methods is that they do 
not need the electromagnetic counterpart to be identified. Classical 
lensing tests can be developed with knowledge of either the source 
redshift or the distance. This is different from the usual cosmo- 
graphic approach proposed for LISA, which attempts to build up 
the Hubble diagram from the measured distances to the binaries. 
Inspiral GWs encode the luminosity distance to a binary but they 
do not encode the source cosmological redshift. To build the Hub- 
ble diagram, the electromagnetic counterpart of the GW emission 
needs to be localised independently to determine the event redshift. 
Lensing methods can overcome such a shortcoming. 

LISA sources are expected to lie at very high redshifts. Even 
if the accuracy of the proposed lensing tests is not competitive with 
results from type la supernovae or the cosmic microwave back- 
ground, the explored redshift range would be quite unique and cos- 
mological parameters might be tested with direct methods out to 
z < 10. The proposed methods are still based on cosmological 
distances but, due to the distinct redshift range explored, are or- 
thogonal to other tests at lower z. 



This circumstance could make very interesting their use in 
combination with other techniques. Let us just consider a still very 
mild constraint on the dark energy equation of state, whose upper 
allowed limit is lowered from to —0.5 in the following example. 
In our optimistic scenario (3 lenses, SD^t / D^t = 10%), the ac- 
curacy on the dark matter parameter from the time-delay test would 
improve from SQm — 0.11 to 0.09 for unknown source redshifts 
or from SQm — 0.08 to 0.06 for known z s 's. 

At the very least, even in a pessimistic scenario with just one 
lens the cosmological concordance model might be tested in an in- 
dependent context and the Hubble constant would be directly mea- 
sured up to very high redshifts. The ability to identify the lens is 
crucial to the discussed methods. For high redshift sources with 
a small signal-to noise ratio, the sky location accuracy might be 
quite poor, ~ ldeg 2 . Several hundreds of sources are expected 
to be localised in such error-box. However, to claim lensing two 
above threshold signals have to be observed in the same sky posi- 
tion. By combining the statistical estimated positions for the two 
sources, the final error-box should be reduced by half. Lensing the- 
ory can also contribute to candidate selection. Deflectors should be 
around redshift one, where we should consider a comoving volume 
of ~ 1 — 3 x 10 10 Mpc 3 . Lenses are also expected to be very mas- 
sive, which makes them rare. Considering a numerical density of 
> 10~ 4 Mpc -3 , we expect just a dozen of lens candidates in the 
LISA error box. 

In a more optimistic scenario, the lensed signal might be very 
loud, which would allow a much more precise determination of the 
source position. In a best case scenario, lensing amplification of the 
source could be enough to make the electromagnetic counterpart of 
the source detectable in both images. 

The certified accuracy of cosmological lensing methods with 
LISA can be obtained only with a good knowledge of the galaxy 
formation history and merging rates. This is the essential ingredi- 
ent to estimate the number of multiple lensing events expected to be 
detected by LISA iSereno et al.ll201oT) . Nevertheless, even in pes- 
simistic scenarios cosmological tests based on gravitational lensing 
of GWs by LISA seem to offer a new and independent way to in- 
vestigate the properties of dark matter and dark energy and fix the 
global distance scale. From an alternative point of view, if we are 
very confident on the estimates of cosmological parameters as ob- 
tained from independent tests, lens number counts with LISA might 
be used to constrain galaxy formation and evolution up to z < 3. 
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